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Abstract 

"Decoherence of quantum superpositions through coupling to engineered 
reservoirs" is the topic of a recent article by Myatt et al. [Nature 403 , 269 
(2000)] which has attracted much interest because of its relevance to current 
research in fundamental quantum theory quantum computation, teleporta- 
tion, entanglement and the quantum-classical interface. However, the pre- 
ponderance of theoretical work on decoherence does not consider the effect of 
an external field . Here, we present an analysis of such an effect in the case of 
the random delta-correlated force discussed by Myatt et al. 
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"Decoherence of quantum superpositions through coupling to engineered reservoirs" [1] 
is the topic of a recent article which has attracted much interest because of its relevance 
to current research in fundamental quantum theory, quantum computation, teleportation, 
entanglement and the quantum-classical interface. As Schleich remarks in an accompanying 
"News and Views" article, this is "- a pioneering experiment that engineers decoherence -" 
[2]. However, the preponderance of theoretical work on decoherence [3,4] does not consider 
the effect of an external field . Here, we present an analysis of such an effect in the case of 
the random delta-correlated force discussed in Ref. 1. 

Myatt et al., [1] used a linear Paul trap to confine single Be ions in a harmonic potential 
and then prepared various superposition states. Next, they induced decoherence by coupling 
the single ion to a reservoir which they controlled in various ways. Such a reservoir gives 
rise to an external force f(t) in the equation of motion of the system, in contrast to the 
usual intrinsic fluctuation force F(t) which arises from interaction with an ambient thermal 
dissipative environment [4], which, of course, will always be present, even at T = 0. Thus, 
the question arises as to not only what is the dependence of the characteristic decoherence 
decay time r d on the separation d of the superposition components, the temperature T 
and the dissipative decay rate 7 (all of which come into play when f(t) = 0, the focus of 
most theoretical work) but also what is the dependence on the parameters of the engineered 
reservoirs which give rise to f(t). The existing experiments [1] focused on the dependence 
of r d on d and demonstrated that r d ~ d~ 2 . This is a familiar result predicted by the 
plethora of papers dealing with the f(t) = situation but it does not give information 
on the dependence of r d on the parameters of the externally-superimposed reservoir. More 
details on the experimental results were given by Turchette et al. [5] and these authors 
also reviewed the theory of the damping of a harmonic oscillator in a dissipative reservoir. 
Whereas the latter gives rise to a fluctuation force on the oscillator which is related to 
the dissipation via the fluctuation-dissipation theorem, an externally engineered situation 
requires an additional analysis, as is made clear in [6]. 

Much of the discussion of decoherence [3,4] has been in terms of a particle moving in 
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one dimension that is placed in an initial superposition state (a Schrodinger "cat" state) 
corresponding to two widely separated Gaussian wave packets. The corresponding wave 
function has the form 

^ 0) = (87r^)V4(i I ( ex P {-^4^} + 6XP {"^4^}) ' (1) 

where d is the separation and a 2 is the variance of each packet. It is clear that the spatial 
probability distribution P(x, 0) for this superposition of two states consists of the sum of the 
probability distributions for the individual states plus an interference term. In the absence 
of dissipation [7,8], one proceeds by calculating ip(x,t) from which P(x,t) readily follows. 
However, when dissipation is present it is necessary to use a density matrix approach [4] 
which, when combined with the use of quantum probability functions, led to an expression 
for P(x, t) of the form 

P(x, t) = Pi (:r, t) + P 2 (x, t) + P^x, t) cos[/(t)], (2) 

where Pi and P 2 correspond to the time dependent probability distributions for the separate 
wave packets and the third term is an interference term. The latter is characterized by a 
cosine factor (which varies in time according to a known function fit) [4]) which is multiplied 
by an amplitude factor Pj(x,t), which is found to decay in time. The disappearance of the 
interference term, that is the decoherence, is measured by defining an attenuation coefficient 
a(t), which is the ratio of the factor multiplying the oscillatory term to twice the geometric 
mean of the first two terms, i.e. 

a{t) = m . (3) 

2[P 1 (x,t)P 2 (x,t)} 1/2 

We should mention that, in the literature, one finds various measures of decoherence, based 
on decay of diagonal and off-diagonal density matrix elements or probability distributions 
in phase space, momentum space or coordinate space [9] but we consider the latter to be 
the most desirable because it is closest to experiment. Thus, returning to (3), what we have 
found [4] is that a(t) depends crucially on the spreading of the wave packets corresponding 
to the individual states. 



For f(t) = 0, this spreading arises from the possible intrinsic spreading associated with 
the uncertainty principle, thermal spreading and spreading due to dissipative (7) effects. 
Explicitly, for a free particle described by a single wave packet, the width after a time t is 
w(t), given by [4] 

_ [*(«,), X(t, +t )]' 
4<7 Z 

where a is the initial width and s (t) is the mean square displacement (discussed in more 
detail below). For the attenuation coefficient in the case of a free particle we have the 
formula [4] 

In addition, the characteristic time for decay to occur, say, is defined as usual [4,7,8] as 
the time at which a(t) = exp(— 1). 

We now turn to the case where fit) 7^ and we generalize from the case of a free 
particle to that of an oscillator potential, corresponding to the experiment described in 
[1]. For f{t) 7^ 0, there is an additional spreading of the wave packets, which we will now 
calculate. Afterwards, we will turn to the role it plays in the calculation of a(t). 

Let x(t) be the dynamical variable corresponding to the coordinate of the wave function 
of the superposition state of the oscillator of Myatt et al. [1]. As shown in Ref. [10], in the 
presence of an external force f(t) in addition to the fluctuation force F(t), the steady-state 
motion can be described by means of a generalized quantum Langevin equation 

mx+ f dt'n(t-t')x(t') + Kx = F(t) + f(t), (6) 

where /i(t) is the memory function, K is the oscillator force constant (K = mu^), where u 
is the oscillator frequency, and F(t) is a fluctuating operator force with mean (F(t)} = 0. 
The steady-state solution of (6) can be written as 

x(t) = f dt'G(t - t')[F(t') + f(t')} 

= x s (t) +x d (t), (7) 



where x s (t) is the stationery solution and Xd is due to the driven motion. Also, G(t), the 
Green function, is given by 

1 f°° 

G(t) = — duja(uj + iO + )e- tuJt , (8) 

2,71 J-oo 

with a(z) the familiar response function 



—mz 2 — izn(z) + K 



In addition 



roc 

ji(z) = / dtfi(ty zt 

Jo 

= m 7 (z), (10) 

is the Fourier transform of the memory function and it characterizes the reservoir [4,10]. 
The fact that ct(z) does not depend on f(t) follows simply by taking the Fourier transform 
of (7) which enables the solution to be written in Fourier transform language as 

x(u)=a(u)[F(u) + f(u)], (11) 

where superposed tildes indicate Fourier transforms. 

Because of the linearity of the oscillator, it is clear that the motion of the driven oscillator 
will be a superposition of a driven mean motion and a motion about the mean that is identical 
with the motion about the equilibrium state [6]. The starting-point of our calculation is the 
correlation 

l -(x(t)x(t') + x(t')x(t)) = C(t - t') = C + C d 

= -/ dtolm{a(uj + i0 + )}coth— — cosu(t-t r ) + C d , (12) 
7r Jo 2kT 

where C and Cd are the contribution due to F(t) and f(t), respectively. It follows that the 

mean-square displacement (which characterizes the spreading of the wave packet) is 

s(t) = (l X (t)-x(0)f) = 2{C(0)-C(t)} 

= — / dtulm{a(cu + i0 + )} coth — — (1 — cosut) + Sd, (13) 
7i Jo ' 2kl 



where C(t) is given by (12) and s d is the contribution due to the "driven motion". 

Here, we have used the fact that since < F(t) >= and since there is no correlation 
between F(t) and /(£), it is clear that 

s(t)=s (t)+s d (t), (14) 

where Sq denotes the contributions due to F(t). Since So(t) has been calculated in detail, 
in Ref. [4], which considers entanglement between the system and the environment at the 
initial time t — 0, we will henceforth concentrate on s d . Consider that the external force is 
applied at t — 0. It follows from (7) that 

x d (t) = f dt'G(t-t')f(t'). (15) 
J o 

Since Xd(0) = 0, it follows that 

s d (t) = (x 2 d (t)} 

= f dt' f dt"G(t - t')G(t - t")g(t' - t"), (16) 
Jo Jo 

where 

g(t'-t") = (f(t')f(t")). (17) 

Further progress clearly depends on the nature of f(t) but, keeping in mind the existing 
experiments [1,5], let us consider a random delta-correlated force so that 

g (t>-t")=g6(t'-t"), (18) 

where g is time-independent. Hence, substituting (18) in (16), we obtain 

s d = g f dt'G 2 {t'). (19) 
Jo 

In the case of the oscillator potential of Myatt et al. [1], we find that in the case of Ohmic 
coupling (7(0;) = 7 = constant) 

G(t) = e-^ (20) 
muji 
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where 

uj\ = u 2 - ( 7 /2) 2 . (21) 
Thus, substituting (20) in (19), it follows that 

s d = Am2 9 luj2uJ 2 {(! - e" 7t )2^ 2 - e- 7t (7 2 sin 2 uj Y t + 7^ sin 2^)} • (22) 

In the absence of dissipation (7 — > 0) (which approximates the experiment of Myatt et al.), 
(21) and (22) give 

g f sin2w t\ , . 



2m 2 Uo I 2cj t J 

Again, for 7 — > and T — > (absence of dissipation and for negligibly low temperatures), 
it readily follows that 

where a is the initial width of the individual wave packets. Thus, the dependence on d 2 in 
the numerator always emerges, regardless of the value of s d . We also note the absence of a 
term analogous to the second term in (4), corresponding to the fact that, when f(t) = 0, 
the width of the oscillator wave function is constant in time whereas that of the free particle 
continually increases. 

It is clear from (24) that the relative magnitudes of s d and the initial variance a 2 play a 
crucial role. In particular, 
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a(t) w expj-^ j> ifs d >>a 2 , (25) 



and 



s d d 
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a(t) w exp j-— — I if s d « a 2 . (26) 

Thus, in the former case, the result for a(t) is independent of s d i.e. independent of the 
external force f(t). In the latter case, using (23), we see that 



where 



16a 4 m 2 ^ 2 

T ° = (28) 



For small times (2uj t « 1) after the initial time t — 0, we see that 

r i (2cj t) 2 l 



a(i) = exp < — 



r 6 
gd 



2 



= exp i-»' 3 j- if "° t<<1 ' (29) 

in which case the decay rate of decoherence is independent of u>o, corresponding to free 
particle behaviour. However, when t further increases there is a change in the time behaviour 
until at the end of the first cycle at 2u t = 2n, we see from (27) that 

a(t) = expj- — J. (30) 

In fact, as we go into the next and subsequent cycles, the sin(2u;o^)/2u;o^ term becomes more 
and more negligible so that (30) becomes more and more accurate as we go beyond the first 
cycle. 

It should also be noted that (18) also corresponds to a white-noise spectrum. However, 
it is very different in nature than the white-noise spectrum associated with the fluctuation 
force F(t). A random c-number field feeds energy into the quantum particle (and, in fact, 
for a particle with negligibly weak coupling to a heat bath and for either a zero or oscillator 
potential, it may be shown that the energy of the particle increases linearly in time). On the 
other hand, in the case of a fluctuation force, we are necessarily dealing with a heat bath; in 
other words, we have a dynamical system in which the particle also loses energy due to the 
emission of bath excitations. Thus, for example, in the case where the white-noise spectrum 
is associated with an equilibrium temperature [10] 

(F(t')F(t")} = 2m 7 T<5(t' - t"). (31) 
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for the case of constant 7 and in the classical limit. Moreover, the rate of work being done 
by the fluctuation force, Pp say, is given by [11] 

p F = fcry (32) 

Thus, the rate of work being done by the fluctuation force is proportional to the dissipation. 
This is a manifestation of the general principle that, at equilibrium, the energy lost by 
a particle due to dissipation is compensated by the energy received from the fluctuation 
force. Thus, there is a crucial difference between the effects of f(t) and F(t) so that, in 
particular, an external field that has a white noise spectrum can not be approximated by a 
weakly-coupled thermal reservoir and, as a result, one must use the analysis given above. 

Finally, it is clear that in order to explore the larger parameter space (such as dependence 
on T, 7 and various choices of f(t) as well as on the potential), both further experiments 
and theoretical work will be needed. Some recent work has made inroads into this multi- 
dimensional parameter space. First, for the problem considered above, we find that a non- 
random external force does not cause decoherence. Second, in the absence of an external 
field, the Schrodinger cat superposition has been examined for the case of an oscillator 
potential and high temperature [12] and for the case of a free particle subject to the effects 
of the zero-point oscillations of the electromagnetic field [13]. 

We are pleased to thank Professor G. W. Ford for many enlightening discussions. 
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